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Abstract 

The interaction of two-level atoms with a common heat bath leads to an effective interaction 
between the atoms, such that with time the internal degrees of the atoms become correlated or 
even entangled. If part of the atoms remain unobserved this creates additional indirect decoherence 
for the selected atoms, on top of the direct decoherence due to the interaction with the heat bath. 
I show that indirect decoherence can drastically increase and even dominate the decoherence for 
sufficiently large times. I investigate indirect decoherence through thermal black body radiation 
quantitatively for atoms trapped at regular positions in an optical lattice as well as for atoms 
at random positions in a cold gas, and show how indirect coherence can be controlled or even 
suppressed through experimentally accessible parameters. 



1 



With the rise of quantum information processing it has become necessary to understand 
decoherence in true many-particle systems. It has been known for a long time that the 
decoherence rate for a single degree of freedom scales like a power of a certain distance 
between the components of a superposition. The "distance" and its natural scale depend on 
the coupling to the heat bath. For example, if the single degree of freedom couples through 
a spatial coordinate x to the heat bath, the latter selects eigenstates of x as "pointer-basis" 
Q], and the relevant "distance" is measured in configuration space, with a microscopic length 
scale such as the thermal de Broglie length as natural unit. Decoherence therefore becomes 
extremely fast for mesoscopic or even macroscopic distances, and this is considered one of 
the main reasons why the everyday world around us behaves classically. Other couplings 
lead to different power laws and different natural microscopic units [2]. Decoherence pro- 
cesses for single degree of freedom systems are nowadays routinely resolved experimentally 
for microscopic distances between the superposed components, and have been so far in good 
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agreement with the theoretical predictions Pl y, 1^ , lO, 13, ISi • However, decoherence measure- 
ments on true many particle systems are only now becoming available j^, and there is a need 
for detailed theoretical predictions, in order to verify the validity of quantum mechanics in 
an entirely new regime, namely one where the joined states of many particles are coherently 
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Theoretical progress was achieved recently with the derivation of a "decoherence metric" 
jl^j ]. which measures the distance between the components of a quantum superposition of 
arbitrarily many qubits with degenerate energy levels, and determines directly the time de- 
pendent decoherence. It turned out that for sufficiently far separated qubits with degenerate 
energy levels the time dependent decoherence boils down to just single qubit decoherence 
multiplied with the standard Hamming distance between the superposed quantum code 
words. For smaller qubit separations (the relevant length scale is the inverse of the wave 
length of the UV cut-off of the bath modes), interference effects start to play a role and one 
sees strong deviations from the simple scaling with the standard Hamming distance. Never- 
theless, the notion of a distance (more precisely: a pseudo-metric, in the strict mathematical 
sense), can be maintained through the introduction of a metric tensor determined by the 
heat bath, whose off-diagonal elements reflect the interference processes. As a consequence, 
2N{N-i)/2 independent decoherences are governed by only ~ A^^ matrix elements of the met- 
ric tensor. The entanglement of a state alone does not determine how fast it decoheres: For 
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example, a GHZ state (|000) + |lll))/v^ has maximum Hamming distance between its two 
components, but if the qubits are sufficiently close the decoherence metric will distinguish 
this state for example from the state (|001) + |110))/-\/2, or all other states which differ from 
GHZ by flipping qubits, whereas all theses states have the same entanglement. 

The metric tensor contains the contribution from the direct decoherence process discussed 
so far, but also an "indirect decoherence": the heat bath generates effective interactions 
between the qubits which can lead to classical correlations or even entanglement between 
them ("reservoir induced entanglement", ^3). "Indirect decoherence" is the additional 
decoherence that is induced if some of the atoms which got correlated or entangled with the 
selected atoms remain unobserved. In this paper I investigate indirect decoherence in more 
detail, and show that even for a rather small number of unobserved atoms (of the order of 
10) and for weak effective coupling between the atoms, indirect decoherence can drastically 
enhance the overall decoherence. The effect should be important if one wants to build a 
quantum memory from trapped atoms in an optical lattice, but when additional atoms get 
trapped in the optical lattice and are not read out. Indeed, until recently Q] it was difficult 
to even control the number of atoms per lattice site. The situation might be even worse 
for quantum information stored in macroscopic gas samples p], where the total number of 
atoms in which the information is stored can only be estimated and one has no control over 
which individual atoms store the quantum information. 



I. THE MODEL 



Let us consider two level atoms at arbitrary but fixed positions Rj {i = 0, . . . N — 1) 
interacting with thermal black body radiation, which forms a common heat bath. All atoms 
are assumed identical with degenerate energy eigenstates | — 1) arid 11), with cx^lil) = ±|±1). 
In dipole coupling approximation, the total hamiltonian reads [l6 | 

N-l 

k k i=0 

where a^i and a^i are Pauli matrices for atom i. The index k stands for wave vector k and 
polarization direction A [kj = 2Ttnj/L with integer rtj, j = x,y,z for periodic boundary 
conditions); a|. (ofc) are the creation (annihilation) operators for mode k with frequency 



iUk = c|k|, polarization vector e^, and electric field amplitude £ = \/hwkl (2^0^) , where c, 



and V are the dielectric constant of the vacuum, speed of hght, and the quantization volume, 
respectively. The coupling constant of atom i to mode k is denoted by g^^^ = —^vS^'^ ■ e^, 
where u*^*^ stands for a unit vector in the direction of the dipole moment of atom i, (— l|d|l) = 
erfw^'^ with electron charge e and dipole length d. The restriction to atoms with degenerate 
energy levels, VLq = 0, leads to a vanishing system hamiltonian, ifgys = |^^o Xlilo^ '^zi = 0. 
The model is a special case of the more general class of models known as pure dephasing 
models, where the system hamiltonian commutes with the interaction hamiltonian, i.e. the 
second term in (Q). These models can be solved exactly for an arbitrary number of atoms 
at arbitrary positions. 



II. DECOHERENCE METRIC 



We are interested in the decoherence process of the n selected atoms (indices 0, . . . , n — 
1) out of the atoms. We therefore have to first trace out the electro-magnetic (e.m.) 
field modes, leaving a density matrix p, and secondly the unobserved atoms n. . .N — 1. 
The resulting reduced density matrix p{t) of the remaining atoms will be expressed in the 
eigenbasis of the axi, axi\^)x = the natural basis (also called pointer basis) for 

studying the decoherence process It has matrix elements Pss'{t) = tr„..jv-iPss'(^), where 
s and s' are subsets of length n of the labels s = (sq, si, . . . sn^i) and s' = (sq, s[, . . . s'j^^i) 
of the quantum states |s) and |s') of all atoms, taken as column vectors, and Si,s'^ = ±1, 
i = 0,l,...,A^ — 1, refer to atom i. We assume that all unobserved atoms are initially in 
the energy eigenstate |1) = (|1)^ + | — l)x)/V^, and that there are no initial correlations 
between the unobserved atoms and the selected atoms. 

The dynamical quantities of interest are the "decoherences" dss>{t), which we define as 
normalized complements of "coherences" (i.e. off-diagonal elements of the reduced density 
matrix of the n selected atoms alone). 







dss'it) ^ 1 - for p,,,(0) ^ . (2) 

iPss' 



it was shown that the behavior of the decoherences is given by dss'{t) ~ ||s — s'||^^^^^ 
with the "decoherence metric" 

||§-§1|M(i) = ^A/(§-§0^M(t)(s-s'), (3) 



4 



where ^ denotes the transpose, and M(t) is a real, symmetric, and non-negative time 
dependent "decoherence metric tensor" (DMT) with matrix elements = 0, . . . ,n — 1, 
j3 = l/ksT is the inverse temperature), 

M,,- (t) = 4^- (t, Ri - R,) + 2$,,- (t, R,, R,) , (4) 
„(*)„(i) , 
fi,{t, R) = V ^-^^^ cos(k ■ R)(l - coscufct) coth (5) 

(0 (i) 

(t , R) = 2 V ^*L^ cos(k ■ R) {uokt - sin w^t) (6) 

k 

N-l 

^ij{t, Hi, Rj) = v^jfc(^, Rj — Hk)'^jk{t, Hj — Rfc) . (7) 

k=n 

The heat -bath itself therefore induces a natural distance ||s — s'||7v/(t) between the n-qubit 
states, which determines directly the time dependent decoherences. The validity of eq.Q is 
limited to \M,j{t)\ < 1 Vi, j. 

The distance | |s — s'| generalizes the well-known Hamming distance D^{s,s'), which 
is defined as the number of bits in which s and s' differ, and which is obtained for M = I. 
This limit is reached for sufficiently large separation of the qubits, and dgg' (t) then goes over 
into D^{s,s') up to a time dependent function describing single qubit decoherence jl j |. 

It is clear from the definition that M(t) is real and symmetric. In appendix I show 
that M(t) is also non-negative and obeys the triangle inequality. However, if a decoherence 



free subspace (DFS) [itIIisI IiqI l2c]. l2l|| exists, the decoherence metric is strictly speaking a 



pseudo-metric, as there can be code words s and s' with s 7^ s' such that ||s — s'||Af(t) = 0. 



III. ONE SELECTED ATOM 

Indirect decoherence, i.e. the part $jj in eq.(jH) is best appreciated for just n = 1 selected 
atom (taken to have index i = 0), and — 1 non-observed atoms. Only one decoherence 
is then relevant, di_i, and M has a single matrix element, Moo(t) = 4/oo(t) + 2$oo(t). We 
express all lengths in terms of the dipole length d, = jR, — Tij\/d, and times in units 
of d/c. We will furthermore assume that all dipoles are oriented in the same direction u. 
The angle between u and the vector — will be called 6ij, such that (y9jj(t,Rj — Rj) 
becomes a function of t, r^j, and 9ij, and depends on i and j only through these variables, 
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in 0, 



(Pij{t, Rj — Rj) = ^{t, Tij, 6ij). As shown in |13|, /oo is given for T = by 

„ 2 /k^ 1 — cos(K.t) — K.tsin(K,t)\ , . 

w = ^ — H — ^ ) • 

where k = k^^xd is a UV cut-off of the heat bath. A necessary condition for the dipole- 
couphng approximation is k <^ 1. Corrections to (jH)) due to finite temperature are of order 
ksT / (huj^ax) with tUmax = c/cmax, and will be neglected in the following. Eq.® implies that 
for a finite UV-cutoff the direct decoherence will remain finite for all times, 

/oo ^ ^ . (9) 

Such a behavior has been termed "incomplete decoherence" Q|. The initial behavior is 
quadratic in t, /oo ^ 7^t^ with 7 = {12'!t)k^ for 7^ <^ 1. 
The function ip{t,r,9) reads 

ip{t,r,e) = ^|^(l + 3cos(2^)) (si((r + t)/€) -Si((r-t)/€)) 

+ ^^^^^^7—^ ((r + t)Si((r + t)K) - (r - t)Si((r - t)^) - 2tSi(/«r)) | . (10) 

Rapidly oscillating terms of the type sin(/tr), cos(Kr), and Krcos(Kr) have been neglected 
here, as their average in the case of a small uncertainty in the atom positions is exponentially 
small: atoms of mass M trapped in the ground states of harmonic oscillators with trapping 
frequency v have a Gaussian distribution of their center of mass with a width 5r ~ ^JhjvM^ 
leading to a suppression of these terms by a factor exp(— (fi;(5r)2/2). A typical experimental 
parameter, v ~ 30kHz j^], leads to brd ~ lOOnm. Optical dipole lengths d ~A and a UV 
cut-off /cmax ~ 1/loA give K ~ 0.1, and Kbr ~ 100, so that these terms can indeed be safely 
neglected. For a cold gas 5r is expected to be even larger. 
In the limit of |r ± ^ 1, V9jj(t, r, ff) approaches 

^{t, r, B) = a^{3 cos^ 6 - l)0(t/r - 1), (11) 

where B(t/r — 1) is the Heaviside function centered on the light cone. We therefore recognize 
(y9 as a phase accumulated due to an effective dipole interaction between the atoms mediated 
through the modes of the electromagnetic field. Note that in this limit indirect decoherence 
becomes basically independent of the cut-off k. 
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A. Optical lattices 



In the following we consider specifically the situation for a 2D square optical lattice with 
lattice constant a (taken in units of the dipole length d as well), with the single selected 
atom at the center of the lattice. The fact that ip{t^ r,6) oc t leads to an unbound quadratic 
growth of $ij(t,Ri,Rj) with t, $ij ~ Nnn{,Oit / a^Y , where Nnn is an effective number of 
nearest neighbors, weighed by the inverse cube of their distance from the selected atom in 
units of the lattice spacing. As a/t^ -C 1 an immediate consequence of eqs. ()lll9p is that for 
large enough times, t > ti with 

h ~ P^-^ , (12) 

indirect decoherence always dominates over direct coherence even for a small number of 
unobserved atoms close to the selected atom. Due to the strong r-dependence of (fTTj) . the 
nearest neighbors and next nearest neighbors give the by far leading contributions to the 
indirect decoherence. Indirect decoherence dominates immediately (i.e. as soon as t ^ a ^ 
1, where the quadratic behavior of $jj ~ A^n„(at/a^)^ in t is valid) over direct decoherence, if 
Nnn{,oit/a?Y > 7^t^, or a < ttc with the critical spacing Oc ~ (127raA^„„)^/^/fi;^/^. If the cut- 
off is of the order /zo'max ~ leV and d = lA, acd is of the order lOOnm. If the cut-off is given 
by the break-down of the dipole approximation (fcmax ~ 271 /d), ac reduces to 0^, ~ 1, i.e. the 
atoms will have to become basically closely packed before indirect decoherence immediately 
dominates over direct decoherence. 

One might object that a small number of two-level atoms with which the selected atoms 
effectively interact cannot constitute a real heat-bath, and should rather lead to repeating 
revival phenomena of the coherences instead of to decoherence. However, it turns out that 
even for a square optical lattice of 3x3 atoms (i.e. 8 unobserved atoms), the revivals are 
hardly visible, and in a square optical lattice of 31x31 atoms all revivals seem to have 
disappeared completely. This is shown in Figure ^ where we see the decoherence as function 
of time for a 2D square optical lattice with lattice constant a = 1000, dipole moments 
perpendicular to the plane of the lattice, and with the selected atom in the center of the 
lattice. Direct decoherence sets in immediately and increases oc t"^ for small times (see 
eq.®), before saturating at aK^/(37r). Due to the weakness of the dipole coupling, indirect 
decoherence becomes appreciable only at much later times for atoms separated thus far. But 
because of the continued quadratic growth of the indirect decoherence it finally destroys all 
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FIG. 1: The decoherence for a single selected atom in the center of a 2D square optical 

lattice of 31x31 atoms with lattice constant a = 1000 (in units of the dipole length) as a 
function of the dimensionless time t for k = 0.01 (black circles), k = 0.1 (blue squares), and 
K = 1 (red diamonds) along with the decoherence metric predictions (continuous lines). The 
first rise corresponds to the contribution of direct decoherence due to the interaction with the 
e.m. modes, the second rise results from indirect decoherence due to the effective interaction with 
the non-selected atoms mediated by the e.m. modes. Superpositions are decohered completely 
when di-i = 1 is reached. 



coherence left by the direct decoherence. The figure also shows that the exact result for di_i 
jlS^ is very well approximated in the entire interesting regime di-i < 1 by the decoherence 
metric prediction, eq.(jni). 

Indirect decoherence has an interesting dependence on the orientation of the dipoles. 
Eq. flllj) shows that $oo oc (3cos^^^ ~ 1)^- Indirect decoherence in a 2D optical lattice can 
therefore be completely suppressed by orienting the atomic dipoles at the magical angle 
9 = arccos(l/-\/3) — 54.7° with respect to the lattice. This angular dependence might serve 
as additional experimental signature of the effect. 

B. Cold gases 

For an atomic gas, the positions Rj of the atoms are not known. We resort to an ensemble 
description, where we average over the positions of the atoms. We assume that all Hj with 
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the exception of Rq = are randomly, independently and evenly distributed with an average 
density (atoms per volume) py- We find 

N-l 

($oo(t,0,0)) = (5^^^^(t,Rfc)) 



k=l 



= 2TTpvd^ J^drr^ d^sm^{at/rY{3cos^ ^ - iyQ{t/r - 1) . (13) 

The lower cut-off is now given by the smallest distance up to which two atoms might approach 
each other, which for a sufficiently dilute gas at low temperature is of the order of the 
scattering length / (taken in units of d as well), if / > 0. This gives 



{%o{t, 0, 0)) ^ j'at' with 7G = « V • ^^'^^ 

Indirect decoherence immediately dominates over direct decoherence for densities pv ^ 
K^P /d^. A UV cut-off /icumax = leV, d = lA, I = lOA gives a critical density of about 10^° 
atoms/m^. Recent experiments on dense Bose-Einstein condensates deal already with similar 
densities The coherence of internal (spin) degrees of freedom of condensed bosons has 
recently been demonstrated [2^, so that indirect decoherence in a cold gas might become 
observable in the near future. For smaller densities indirect decoherence takes over for t > t2 
with 

'^^'^\/^- ^''^ 
The dependence of the indirect decoherence on / gives the interesting perspective to control 
indirect decoherence in a cold gas through a Feshbach resonance, which allows to vary / over 
many orders of magnitude |25| . 



IV. CONCLUSIONS 



I have shown that indirect decoherence due to reservoir induced entanglement between 
degenerate two-level atoms can substantially increase decoherence in an optical lattice or 
a cold atomic gas, compared to the direct decoherence due to the coupling of each atom 
to the e.m. field. For large enough times indirect decoherence in fact always dominates, 
even for only a few unobserved atoms. For sufficiently densely packed atoms the dominance 
of indirect decoherence begins as soon as a light signal has traveled a dipole length. The 
dependence of the indirect decoherence on the orientation of and the distance between the 
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dipoles offers the interesting perspective to control indirect decoherence with easily accessible 
parameters. In a 2D optical lattice, indirect decoherence can be switched off completely by 
orienting all dipoles under a magical angle = arccos(l/\/3) with respect to the lattice, 
and in a dilute, cold atomic gas, one can suppress indirect decoherence to large extent by 
increasing the scattering length / through a Feshbach resonance. 
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V. APPENDIX A: PROOF OF NON NEGATIVITY AND OF THE TRIANGLE 
INEQUALITY FOR THE DECOHERENCE METRIC 

A. Non negativity 

We show separately '^iXifijXj > and — ^ From eqs. (|5l6p we have 

I,] I,] k 



E 



^^i9k 



(»)gik-i?, 



2 

1 - cosujkt (3f\hJk - „ 

7^ coth > 

^k 2 

2 



i,j i,j I I \ i / 

Thus, also -XiMijXj > 0, and M is therefore non-negative. 



B. Triangle inequality 

We define the linear map (p '■ 11^" ~^ I^"; v Mv, where M is a real symmetric, 
non-negative n x n matrix, i.e. v^Mv > for all v G M". We also define the bi-linear 
form (-, ■) : X M" ^ M, (v, w) = V"^Mw, which is not a scalar product, as (v, v) can 
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be zero for v 7^ 0. One can nevertheless prove the Cauchy — Schwartz (C.S.) inequality 
(v,w)^ > (v, v)(w, w): 

Let Vq be the kernel of 0. Thus Mv = = v-^M Vv G Vq. Suppose first that v G Vq 
or w G Vq. Then (v, w) = 0, but also (v, v)(w, w) = 0, as at least one factor is zero. 
Thus the C.S. inequality is trivially fulfilled. Now suppose that v ^ Vq and w ^ Vq. Define 
V = (w, w)v — (w, v)w. We have 

< (v, v) = (w, w)[(w, w)(v, v) — (w, v)^] . (16) 

It is easily seen that (w, w) 7^ if w ^ Vq: Decompose M = Mq + M+ with Mq = PqMPq, 
M+ = (1 — Po)M(l — Pq), where Pq is the projector onto Vq. M+ is the positive part of the 
map, i.e. w^M+w > Vw 7^ 0. Thus, from (w,w) = follows w = or M_|_ = 0. In both 
cases w G Vq. Thus, for w ^ Vq we have (w, w) 7^ 0, and as M > 0, this means (w, w) > 0. 
Therefore eq. ()16p immediately gives the C.S. inequality. 

The proof of the triangle inequality ||s — s"||Af < ||s — s'||a/ + ||s' — s"||j\/ then proceeds 
in the usual fashion. One defines the norm ||x||jvf = a/ (x, x), and the C.S. inequality gives 
||x + y||l,, = (x,x) + (y,y) + 2(x,y) < (x, x) + (y, y) + 2|(x, y)| < | |x| + | |y| |i/ + 
2||x||a/ ||y||j\/ = (||x||j\/ + llylU/)^- The triangle inequality for the distance (jH)) follows from 
here by setting x = s — s', y = s' — s". 
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